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Abstract 

The recently discovered coexistence of incommensurate antiferromag- 
netic neutron scattering peaks and commensurate resonance in underdoped 
YBa2Cu306+z is calling for an explanation. Within the t-J model, the dop- 
ing and energy dependence of the spin dynamics of the underdoped bilayer 
cuprates in the normal state is studied based on the fermion-spin theory 
by considering the bilayer interactions. Incommensurate peaks are found 
at [(1 ± 5)tt,tt] and [vr, (1 ± 5)n] at low energies with 5 initially increasing 
with doping at low dopings and then saturating at higher dopings. These 
incommensurate peaks are suppressed, and the parameter 5 is reduced with 
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increasing energy. Eventually it converges to the [ir, ir] resonance peak. Thus 
the recently observed coexistence is interpreted in terms of bilayer interac- 
tions. 
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The interplay between antiferromagnetism (AF) and superconductivity (SC) in high T c 
cuprates is well-established by now, |I[ but its full understanding is still a challenging issue. 
Experimentally the inelastic neutron scattering (INS) can provide rather detailed informa- 
tion on the spin dynamics of doped single layer and bilayer cuprates. PHlOl An important 



issue is whether the behavior of AF fluctuations in these compounds is universal or not. 
A distinct feature of single layer L^-xSr^CuC^ (LSCO) is the presence of incommensurate 
antiferromagnetic (ICAF) peaks at low energy INS, i.e., the AF scattering peaks are shifted 
from the AF wave vector [71,71] to four points [tt(1 ±5), tt] and [tt, (1±5)tt] (in units of inverse 
lattice constant) with 5 as the incommensurability (IC) parameter, which depends on doping 
concentration, but not on energy. Moreover, ICAF is observed both above and below T c 
in the entire range of doping, from underdoped to overdoped samples. |]-[| In contrast, a 
sharp resonance peak (around 41 meV) is observed in optimally doped bilayer YBa 2 Cu 3 6+a; 
(YBCO) at the commensurate AF wave vector [tt,tt] in the SC state. 0|| Such a resonance 
has also been observed in underdoped YBCO samples with resonance energy scaling down 
with the SC T c , being present both below and above T c . |J Recently, this resonance peak has 
been observed in another class of bilayer SC cuprates Bi 2 Sr2CaCu 2 08 + 5 (BSCO). 0] Such 
a peak has, however, never been observed in LSCO. A very important new development 
is the observation of ICAF in underdoped YBCO in both SC and normal states, with INS 
pattern and doping dependence being very similar (linear in doping for low dopings) to that 
of LSCO. HH However, the IC peak position is energy dependent in underdoped YBCO. 
A challenging issue for theory is to explain the coexistence of this energy-dependent ICAF 
scattering and commensurate resonance in bilayer cuprates. 

Theoretically the ICAF has been interpreted, among others, in terms of Fermi surface 
nesting [pTT] , [T2|] or stripe formation. [Tj| The energy dependence of IC parameter S on energy 
for underdoped YBCO makes the stripe interpretation rather difficult to accept. On the 
other hand, the commensurate resonance peak has been interpreted as due to spin-1 collec- 
tive (particle-hole) excitations, [P^ , P ,I5|| or particle-particle excitations, JTB| or interlayer 



tunneling. |L7| These theoretical treatments are mostly addressing the SC state, and heavily 



rely on adjusting band structure parameters, like the next nearest neighbor hopping t', etc. 
To the best of our knowledge, the ICAF and commensurate resonance peak in underdoped 
bilayer cuprates have not yet been treated from a unified point of view for the normal state. 
No explicit predictions on doping and energy dependence of the ICAF peaks have been made 
so far. 

In our earlier work using the fermion-spin theory, [|19|] the dynamical spin structure 



factor (DSSF) has been calculated for LSCO within the single layer t-J model 20|| , and the 



obtained doping dependence of the IC parameter 5 is consistent with experiments. [[[H3|] In 
this paper we show explicitly if the bilayer interactions are included, one can reproduce all 
main features in the normal state observed experimentally on YBCO, |||| including the 
doping dependence of ICAF at low energies and [ir, it] resonance at relatively high energy. 
The bilayer band splitting in BSCO has been observed in the angle-resolved photoemission 



spectroscopy in both normal and superconducting states. |L8| The convergence of ICAF 
peaks at lower energies to commensurate resonance peak at higher energy is rather similar 
to the scenario argued in Ref. |TUj for the SC state, and the DSSF we derive from the simple 



t-J model (without additional terms and adjustable parameters) demonstrates explicitly this 
convergence. 

The t-J model in bilayer structures is expressed as, 

H = —t^2i CaiaCai+ija — t± ^2,{C\ i(T C 2ia + ll.C.) — \L ^ C\ ia C a i a 
air/a ia aia 

+ J X/ ' ^ai+r) + J± ' ^ 2i ' ("0 

air) i 

where fj = ±x, ±y, a = 1,2 is plane indices, and S a i = C ai aC a i/2 are spin operators 
with a = (a x , <jy, <j z ) as Pauli matrices. The t-J Hamiltonian is supplemented by the single 
occupancy local constraint Cl ia C aia < 1. This local constraint can be treated properly 



in analytical form within the fermion-spin theory [19| based on the slave particle approach, 

C a i-\ = hliS ai , C ai i = h^S^i, (2) 
where the spinless fermion operator h ai keeps track of the charge (ho Ion), while the pseu- 



dospin operator S a i keeps track of the spin (spinon), and the low-energy Hamiltonian of the 
bilayer t-J model (1) can be rewritten in the fermion-spin representation as, 

H = h a i+-f)ha,i(SaiS ai+i ) + S ai S^ i+fl ) + t± '^2(h\ i h 2 i + h2 i h\i){Si i S 2 i + S^S^) 

air) i 

+ fi ^2 h ai h a i + J e fj ^ S a i ' Sai+f) + J-LcE X] ^ u ' ^2j; (3) 

ai air) i 

with J eff = J[(l — p) 2 — <fi 2 ] and Jj_ c fj = J[(l — p) 2 — where p is the hole doping 
concentration, the holon in-plane and bilayer hopping parameters <fi = (h ai h a i + fj) and <p±_ = 
(/ijj/i 2i ), and S£i (Sai) as the pseudospin raising (lowering) operators. In the bilayer system, 
because of the two coupled CuC>2 planes, the energy spectrum has two branches. In this 
case, the one-particle spinon and holon Green's functions are matrices, and are expressed 
as, 

D(i -j,r- t') = D L {% — j,r — t') + r x D T (i -j,r- r'), 
9(i ~ J, T ~ r ') = 9L(i ~j,T- t 1 ) + r x g T (i -j,r- r'), (4) 

respectively, where the longitudinal and transverse parts are defined as, 



D L (i 


-3, 


T 


~r>) 


— ^TS ai {r)S a j\ 


r')>, 


9h{i 


~ 3i 


T 


-r) 


= -(T T h ai (r)hl J ( 


r')), 


D T {i 


— 3i 


T 


-r) 


= —(T T Sai(T)S a ,j 


(r')) 


9r(i 


— 3i 


T 




= -{T r h a i(T)h ] alj 


(r)} 



(5) 

with a a', while r x is the Pauli matrix in the pseudospin space of the layer index. Within 
this framework, the spin fluctuations only couple to spinons, but the strong correlation 
between holons and spinons is included self-consistently through the holon's parameters 
entering the spinon propagator. Therefore both spinons and holons are involved in the spin 
dynamics. The universal behavior of the integrated spin response and ICAF in underdoped 
single layer cuprates have been discussed within the fermion-spin theory by considering 
spinon fluctuations around the mean-field (MF) solution, |2(J where the spinon part is treated 



by the loop expansion to the second order. Following the previous discussions for the single 
layer case, DSSF of bilayer cuprates is obtained explicitly as, 



5(k, uj) = -2[1 + n B (uj)\ [2ImD L (k, u) + 2ImD T (k, u)] 

4[l + n B M](i?i 1) ) 2 ImSg(k,c) 



- (J^y - 5«IteE«(k,a;)]* + [£«ImEg(k, c)] 2 ' 
where the full spinon Green's function, 



(6) 



D-\k, uo) = D^-\k, uo) - E«(k, u), (7) 
with the longitudinal and transverse MF spinon Green's functions, 

4 o) (k,c)=i/2E^ ) /k 2 -(-r ) ) 2 ], 

V 

4 o) (k,c) = i/2E(-i) y+i 5rv^ - («>pn (8) 

V 

respectively, where v = 1,2, and 

ImS^(k,cj) = ImSP(k,cj) + ImS^(k,w), 

ReE^(k, uo) = ReE^k, w) + ReE^ s) (k, u), (9) 

while ImE^(k, u) (ImE^(k, a;)) and ReE^(k, a;) (ReE^(k, ui)) are the imaginary and real 
parts of the second order longitudinal (transverse) spinon self-energy, respectively, obtained 
from the holon bubble as, 

si s) (k,cu) = (W) 2 E E iw(k,p, P >), 

pp' vv'v" 

E?(k,u;) = (l/AQ'E E (-^^^(k.p.p'.w), (10) 
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„ r, 

^fc+P 



IW(k,p,p',c<;) = (Zt[ lp , +p+k + lk - p/ ] +t ± [(-lf + "" + (-1)"+"" 



2 sr } 



1 « ("") 



I _ F£U(k,p,pQ i^Uk, P , P Q 

I , , tM , ,K) , , , t("0 , , >") I ' 



where 7k = (1/Z) J2fj elk fl , Z is the coordination number, 
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BP =B k - J ±cS [ X ± + 2xl(-l)1[e± + (-1)1 
B k = \[{2e X z + Xhk ~ {e X + 2 X % A = 2ZJ cS , 
e = 1 + 2*0/ J eff , e ± = 1 + 4t ± </> ± / J ±eff , 

*£U(k, P, P') = n F $£,)[l - n F (&)] - n B (J$)[n F ($) - n F (^ p/ )}, 

F$A± P, P') = - n F {&)\ + [1 + M^OlM^ ) - M#20], (12) 

n F(£k) anc ^ ^st^fc ) are t ne fermion and boson distribution functions, respectively, and 
the MF holon and spinon excitations, 

& ] = 2Ztxik + V + 2 X ±t ± (-iy +1 , 

(^r=coi+Ai(-ir\ (is) 

with tu 2 = Aa 2 k + A 2lk + A 3 ,A 2 = X llk + X 2 , 
A 1 = aeX 2 ( X /2 + e X z ), 

A 2 = eA 2 [(l - Z)a(e X /2 + X z )/Z - a(C z + C/2) - (1 - a)/(2Z)} 

-aXJ LeS [e(C z ± + xl) + e ± (C± + e X± )/2], 
A 3 = X 2 [a(C z + e 2 C/2) + (1 - a)(l + e 2 )/(4Z) - ae( X /2 + e X 2 )/Z] 

+aAJ ±cfr [e e± C ± + 2CI] + J 2 ±cS (e 2 ± + l)/4, 
Xi = a\J ±cfi [(e ± x + ex±)/2 + ee ± (xl + X z )], 
X 2 = -aAJ ±cff [ee ±X /2 + e L ( X z + C z ± ) + eCj2] - e ± Jl eS /2, 



(14) 



the spinon correlation functions x = (SaiS ai+fl ), X z = (S z ai S z ai+f ^ x± = {SuS^), 
Xl = C = (l/Z^E^iS^S-^,), and C* = (l/Z^^iS^S*.^,), C ± = 

(1/Z) J2fj{S 2 iSi i+fj ), and C z ± = (l/Z) Y^i^u^i+fj) • In order to satisfy the sum rule for the 
correlation function (S^S^) = 1/2 in the absence of AF long range order (AFLRO), a de- 
coupling parameter a has been introduced in the MF calculation, which can be regarded as 
the vertex correction. |^T| All these parameters have been determined self-consistently, as 
done in the single layer case. PU 



At half-filling, the t-J model is reduced to the Heisenberg AF model, and the AFLRO 
gives rise to a commensurate peak at [1/2, 1/2] (hereafter we use the units of [2tt, 2tt]). In 
Fig. 1, we plot DSSF S(k, oS) in the (k x , k y ) plane at doping p = 0.06, temperature T = 0.1J 
and energy u = 0.35 J for t/J = 2.5, t±/t = 0.25, and J±/J = 0.25, which shows that a 
commensurate-IC transition in the spin fluctuation pattern occurs with doping. At low 
energies and lower temperatures, the commensurate peak close to half-filling is split into 
four peaks at [(1 ± S)/2, 1/2] and [1/2, (1 ± 5)/2}. The calculated DSSF S(k,co) has been 
used to extract the doping dependence of the IC parameter S(p), defined as the deviation of 
the peak position from the AF wave vector [1/2, 1/2], and the result is shown in Fig. 2 in 
comparison with the experimental data taken on YBCO (inset). S(p) increases initially 
with the hole concentration in the low doping regime, but it saturates quickly at higher 
dopings, in semi-quantitative agreement with the experimental data. || Apparently, there 
is a substantial difference between theory and experiment, namely the saturation occurs at 
p = 0.10 in experiment, while the calculation anticipates it already at p « 0.05. However, 
upon a closer examination one sees immediately that the main difference is due to the 
appearance of ICAF at too low dopings in the theoretical consideration. The actual range 
of rapid growth of IC parameter 5(p) with doping p (around 4 ~ 5%) is very similar in 
theory and experiment. 

For considering the resonance at relatively high energy we have made a series of scans 
for S(k,u) at different energies, and the result for doping p = 0.06, t/J = 2.5, t±/t = 0.25, 
J±/ J = 0.25 at T = 0.1J and u = 0.5J is shown in Fig. 3. Comparing it with Fig. 1 for the 
same set of parameters except for u = 0.35J, we see that IC peaks are energy dependent, 
i.e., although these magnetic scattering peaks retain the IC pattern at [(1 ± S)/2, 1/2] and 
[1/2,(1 ± S)/2] in low energies, the positions of IC peaks move towards [1/2,1/2] with 
increasing energy, and then the [1/2, 1/2] resonance peak appears at relatively high energy 
(u r = 0.5 J). To show this point clearly, we plot the evolution of the magnetic scattering 
peaks with energy in Fig. 4. For comparison, the experimental result [|TIJ] of YE^CuaOg+x 
with x = 0.85 (p ~ 0.14) for the SC state is shown in the same figure. A similar experimental 
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result || has also been obtained for YE^CusOe+z with x = 0.7 (p ~ 0.12). Although 
these experimental data were obtained below T c , they also hold for the normal state in the 
underdoped regime x < 0.7 (p < 0.12). || The anticipated position oj r = 0.5 J ~ 50 mev 
f22|| is not too far from the peak « 30mev ~ 37mev observed in underdoped YBCO. |9| 
Moreover, the resonance energy uo r is proportional to p at small dopings. We have also made 
a series of scans for S(h, u) at different temperatures, and found both IC peaks and resonance 
peak are broadened and suppressed with increasing temperature, and tend to vanish at high 
temperatures. This reflects that the spin excitations are rather sharp in momentum space 
at low temperatures, compared with the linewidth, and the inverse lifetime increases with 
increasing temperature. Our result is in qualitative agreement with experiments. 

Now we turn to discuss the integrated spin response, which is manifested by the inte- 
grated dynamical spin susceptibility, and can be expressed as, 

I(u;,T) = (l/AO£x"(k,u;), (15) 

k 

where the dynamical spin susceptibility is related to DSSF by the fluctuation-dissipation 
theorem as, x"(k, u;) = (1 — e _/3tJ )S'(k, uj). The results of I(u,T) at doping p = 0.06 in t/J = 
2.5, tjt = 0.25, and J ± /J = 0.25 with T = 0.1J (solid line) and T = 0.2J (dashed line) are 
plotted in Fig. 5 in comparison with the experimental data |23|] taken from YE^CuaOe+z 
(inset), where the dotted line is the function ~ arctan[aiu;/T-|-a3(u;/T) 3 ] with a\ = 6.6, and 
03 = 3.9. These results show that I(u,T) is almost constant for uj/T > 1 and then begin 
to decrease with decreasing uo/T for to/T < 1. It is quite remarkable that the integrated 
susceptibility in the bilayer cuprates shows the same universal behavior as in the case of the 
single layer cuprates, 0] and is scaled approximately as T) oc arctan[aio;/T-|-a3(c<j/T) 3 ]. 
This result is consistent with experiments. ||23[ 

The DSSF in Eq. (3) has a well-defined resonance character, where S*(k, u) exhibits 
peaks when the incoming neutron energy uj is equal to the renormalized spin excitation 
El = (4 1} ) 2 + B^Re^ T (k,E k ), i.e., W{KM = - (J^f - 5«ReSg (k c , u)f = 
{uj 2 — El c ) 2 ~ for certain critical wave vectors k c . The height of these peaks is determined 
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by the imaginary part of the spinon self-energy l/ImE^(k c , uj). This renormalized spin 
excitation is doping and energy dependent. Since ReE^(k, uj) = ReE^(k, u) + ReH^ (k, uj) 
with ReE^ s) (k,a;) < and ReE^fou;) > 0, there is a competition between ReE^ (k, uj) 
and ReEj/(k, u;), which comes entirely from the bilayer band splitting. [18|] At low energies 
the main contribution to ReE^(k, uj) comes from ReE^(k,w), then ICAF emerges, where 
the essential physics is almost the same as in single layer cuprates, and detailed explanations 
have been given in Ref. ||2U|| . Near half-filling, the spin excitations are centered around the 
AF wave vector [1/2, 1/2], so the commensurate AF peak appears there. Upon doping, 
the holes disturb the AF background. Within the fermion-spin framework, as a result of 
self-consistent motion of holons and spinons, ICAF is developed beyond certain critical 
doping, which means, the low-energy spin excitations drift away from the AF wave vector, 
or the zero of W(k$,uj) is shifted from [1/2,1/2] to k<5, where the physics is dominated 
by the spinon self-energy ReE^(k, a;) renormalization due to holons. In this sense, the 
mobile holes are the key factor leading to ICAF. However, ReE^(k, uj) cancels out most 
contributions from ReE^ (k, uj) at relatively high energy, then the anomalous [1/2,1/2] 
resonance reappears. Therefore the bilayer band splitting plays a crucial role in giving rise 
to the resonance. What we calculate is the acoustic spin excitation with modulations in 
the c-direction oc sin 2 (ttzc u L), where zq u is the distance between two nearest Cu layers, L 
the c-axis coordinate in the reciprocal space. This reflects the antiferromagnetic coupling 
between layers, and it is fully confirmed by experiments. fSHlOl 



In conclusion we have shown that if the strong spinon-holon interaction and bilayer 
interactions are taken into account, the t-J model per se can correctly reproduce all main 
features of INS experiments in the normal state in underdoped bilayer cuprates, including the 
doping and energy dependence of ICAF at low energies and [1/2, 1/2] resonance at relatively 
high energy. In fact the ICAF peaks converge to the commensurate resonance, as the energy 
is increased. In our opinion, the difference of AF fluctuation behavior between LSCO and 
YBCO (BSCO) is not due to the presence/absence of stripes, but rather because of the 
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single/double layer structure. Of course, this has to be checked by further experiments. It 
is possible that at some particular energy, a strong commensurate resonance peak coexists 
with the weaker IC features as shown in Fig. 3. 

After submitting this paper, we became aware of the recent INS measurements |24j 



providing evidence for a sharp commensurate resonance peak below T c in the single layer 
cuprate Tl2Ba2Cu6+<5 near optimal doping. However, above T c , the experimental scans show 
a featureless background that gradually decreases in an energy- and momentum-independent 
fashion as the temperature is lowered. The INS in the SC state has not been considered 
so far within the fermion-spin approach, and we need to extend our studies for both single 
layer p0| and bilayer cases to the SC state, where the holon Cooper pairs are formed, and 



the spinon self-energy is originating from both normal and anomalous holon bubbles. Hence 
the renormalized spin excitation in the SC state is very much different from that in the 
normal state, and it may be related to the magnetic peaks detected in the SC state. These 
and other related issues are under investigation now. On the other hand, we emphasize 
that although the simple t-J model can not be regarded as a comprehensive model for the 
quantitative comparison with the doped cuprates, our present results for the normal state 
are in semi-quantitative agreement with the major experimental observations in the normal 
state of the underdoped bilayer cuprates. [Pl,|l0l,|23l . 



ACKNOWLEDGMENTS 

This work was supported by the National Natural Science Foundation of China under 
Grant Nos. 10074007, 10125415, and 90103024. 



11 



REFERENCES 



[1] For reviews, see M.A.Kastner, R.J. Birgeneau, G. Shiran, Y. Endoh, Rev. Mod. Phys. 
70, 897 (1998); A.P. Kampf, Phys. Rep. 249, 219 (1994). 

[2] K. Yamada, C.H. Lee, K. Kurahashi, J. Wada, S. Wakimoto, S. Ueki, H. Kimura, Y. 
Endoh, S. Hosoya, and G. Shirane, Phys. Rev. B 57, 6165 (1998), and references therein. 

[3] G. Aeppli, T.E. Mason, S.M. Hayden, H.A. Mook, and J. Kulda, Science 278, 1432 
(1997). 

[4] J. Rossat-Mignod, L.P. Regnault, P. Bourges, P. Burlet, J. Bossy, J.Y. Henri, and G. 
Lapertot, Physica C 185-189, 86 (1991). 

[5] H.A. Mook, M. Yethiraj, G. Aeppli, T.E. Mason, and T. Armstrong, Phys. Rev. Lett. 
70, 3490 (1993); H.F. Fong, B. Keimer, PW. Anderson, D. Reznik, F. Dogan, and LA. 
Aksay, Phys. Rev. Lett. 75, 316 (1995); P. Bourges, L.P. Regnault, Y. Sidis, and C. 
Vettier, Phys. Rev. B 53, 876 (1996). 

[6] P. Dai, M. Yethiraj, H.A. Mook, T.B. Lindemer, and F. Dogan, Phys. Rev. Lett. 77, 
5425 (1996); H.F. Fong, B. Keimer, F. Dogan, and LA. Aksay, Phys. Rev. Lett. 78, 713 
(1997); P. Bourges, L.P. Regnault, Y. Sidis, J. Bossy, P. Burlet, C. Vettier, J.Y. Henry, 
and M. Couach, Europhys. Lett. 38, 313 (1997). 

[7] H.F. Fong, P. Bourges, Y. Sidis, L.P. Regnault, A. Ivanov, G.D. Gu, N. Koshizuka, and 
B. Keimer, Nature 398, 588 (1999); H. He, Y. Sidis, P. Bourges, G.D. Gu, A. Ivanov, 
N. Koshizuka, B. Liang, C.T. Lin, L.P. Regnault, E. Schoenherr, and B. Keimer, Phys. 
Rev. Lett. 86, 1610 (2001). 

[8] H.A. Mook, P. Dai, S.M. Hayden, G. Aeppli, T.J. Perring, and F. Dogan, Nature 395, 
580 (1998); M. Arai, T. Nishijima, Y. Endoh, T. Egami, S. Tajima, K. Tomimoto, Y. 
Shiohara, M. Takahashi, A. Garret, and S.M. Bennington, Phys. Rev. Lett. 83, 608 
(1999); P. Dai, H.A. Mook, S.M. Hayden, G. Aeppli, T.J. Perring, R.D. Hunt, and F. 

12 



Dogan, Science 284, 1344 (1999). 

[9] P. Dai, H.A. Mook, R.D. Hunt, and F. Dogan, Phys. Rev. B63, 54525 (2001), and 
references therein. 

[10] P. Bourges, Y. Sidis, H.F. Fong, L.P. Regnault, J. Bossy, A. Ivanov, and B. Keimer, 
Science 288, 1234 (2000); P. Bourges, B. Keimer, L.P. Regnault, and Y. Sidis, |conc^ 
mat/ 0006085 . 



[11] See, e.g., N. Bulut, D. Hone, D.J. Scalapino, and N.E. Bickers, Phys. Rev. Lett 64, 
2723 (1990); Q. Si, Y. Zha, K. Levin, and J.P Lu, Phys. Rev. B 47, 9055 (1993); T. 
Tanamoto, H. Kohno, and H. Fukuyama, J. Phys. Soc. Jpn. 63, 2739 (1994). 

[12] See, e.g., G. Blumberg, B.P Stojkovic, and M.V. Klein, Phys. Rev. B 52, 15741 (1995); 
D.Z. Liu, Y. Zha, and K. Levin, Phys. Rev. Lett. 75, 4130 (1995) and more references 
in Ref . [f§ . 



[13] See, e.g., J. Zaanen and O. Gunnarsson, Phys. Rev. B 40, 7391 (1989); D. Poilblanc 
and T.M. Rice, Phys. Rev. B 39, 9749 (1989). 

[14] D.K. Morr and D. Pines, Phys. Rev. Lett. 81, 1086 (1998); F. Onufrieva and P. Pfeuty, 
|cond-mat/99030971 ; A. Abanov and A.V. Chubukov, Phys. Rev. Lett. 83, 1652 (1999). 



JY. Kao, Q. Si, and K. Levin,, Phys. Rev. B 61, 11898 (2000); M.R. Norman, Phys. 
Rev. B 61, 14751 (2000). 



[15] J. Brinckmann and P.A. Lee, Phys. Rev. Lett. 82, 2915 (1999); ibid, |cond-mat/0107138 . 



[16] E. Dernier and S.C. Zhang, Phys. Rev. Lett. 75, 4126 (1995); E. Dernier, H. Kohno, 
and S.C. Zhang, Phys. Rev. B 58, 5719 (1998). 

[17] L. Yin, S. Chakravarty, and PW. Anderson, Phys. Rev. Lett. 78, 3559 (1997). 

[18] D.L. Feng, N.P. Armitage, D.H. Lu, A. Damascelli, J.P. Hu, P. Bogdanov, A. Lanzara, 
F. Ronning, K.M. Shen, H. Eisaki, C. Kim, Z.X. Shen, J.-i. Shimoyama, and K. Kishio, 



13 



Phys. Rev. Lett. 86, 5550 (2001). 

[19] Shiping Feng, Z.B. Su, and L. Yu, Phys. Rev. B 49, 2368 (1994); Mod. Phys. Lett. B7, 
1013 (1993). 

[20] Feng Yuan, Shiping Feng, Zhao-Bin Su, and Lu Yu, Phys. Rev. B 64, 224505 (2001); 
Shiping Feng and Zhongbing Huang, Phys. Rev. B 57, 10328 (1998). 

[21] Shiping Feng and Yun Song, Phys. Rev. B 55, 642 (1997); J. Kondo and K. Yamaji, 
Prog. Theor. Phys. 47, 807 (1972). 

[22] S. Shamoto, M. Sato, J.M. Tranquada, B.J. Sternlib, and G. Shirane, Phys. Rev. B 48, 
13817 (1993). 

[23] R.J. Birgeneau, R.W. Erwin, P.G. Gehring, B. Keimer, M.A. Kastner, M. Sato, S. 
Shamoto, G. Shirane, and J.M. Tranquada, Z. Phys. B87, 15 (1992); B.J. Sternlieb, G. 
Shirane, J.M. Tranquada, M. Sato, and S. Shamoto, Phys. Rev. B 47, 5320 (1993). 

[24] H. He, P. Bourges, Y. Sidis, C. Ulrich, L.P. Regnault, S. Pailhes, N.S. Berzigiarova, 
N.N. Kolesnikov, and B. Keimer, Science 295, 1045 (2002). 



14 



FIGURES 

FIG. 1. The dynamical spin structure factor in the (k x , k y ) plane at doping p = 0.06, tempera- 
ture T = 0.1J and energy u> = 0.35 J for t/J = 2.5, t±/t = 0.25, and J±/J = 0.25. A quasiparticle 
damping T = 0.01J has been used in all results presented. 

FIG. 2. The doping dependence of the incommensurability 5(p) of the antiferromagnetic fluc- 
tuations. Inset: the experimental results on YBCO taken from Ref. ||. 

FIG. 3. The dynamical spin structure factor in the (k x ,k y ) plane at p = 0.06 for t/J = 2.5, 
t±/t = 0.25, Jl/J = 0.25 and u = 0.5J at T = 0.1J. 

FIG. 4. The energy dependence of the position of the incommensurate peaks at p = 0.06 and 
T = 0.1J for t/J = 2.5, t±/t = 0.25, and J±/J = 0.25 (left ordinates) vs the experimental results 
on YBa2Cu306.85 in the superconducting state taken from Ref. [|l0| (right ordinates). 

FIG. 5. The integrated susceptibility at p = 0.06 for t/J = 2.5, t±/t = 0.25, and J L /J = 0.25 
in T = 0.1J (solid line) and T = 0.2J (dashed line). The dotted line is the function 
&iarctan[aiLj/T + ^(uj/T) 5 ] with a\ = 6.6 and a% = 3.9. Inset: the experimental result on 
YBa 2 Cu 3 07_ z taken from Ref. ||. 
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